Time scale of forerunners in quantum tunneling 
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The forerunners preceding the main tunneling signal of the wave created by a source with a 
sharp onset or by a quantum shutter, have been generally associated with over-the-barrier (non- 
tunneling) components. We demonstrate that, while this association is true for distances which 
are larger than the penetration lenght, for smaller distances the forerunner is dominated by under- 
the-barrier components. We find that its characteristic arrival time is inversely proportional to the 
difference between the barrier energy and the incidence energy, a tunneling time scale different from 
both the phase time and the Biittiker-Landauer (BL) time. 

PACS numbers: 03.65.Xp, 03.65.Ta, 03.65.-w 



I. INTRODUCTION 

Tunneling is one of the paradigms of quantum the- 
ory. Textbooks only discuss its stationary aspects, em- 
phasizing the non-vanishing transmission probability for 
monochromatic wavefunctions with energies below the 
barrier maximum. In recent times, since the seminal 
work of Biittiker and Landauer 0, its time dependent 
aspects have been also investigated to find characteristic 
time scales that summarize, together with the transmit- 
tance, the system behavior. Much of this work has been 
controversial, as several authors have proposed and de- 
fended different "tunneling times". In fact each has its 
own virtues, weaknesses, and range of applicability. For 
a recent and rather extensive multi-author review see |2| ; 
for previous reviews see [|[ ||, |(|. An analysis of the 
involved non-commuting observables (the projectors that 
determine the final transmission and the probability to 
find the particle in the barrier region) shows that, from 
a fundamental perspective, there is no unique tunneling 
time, because several quantizations are possible due to 
different operator orderings and defining criteria W . 

This does not mean though that the timing question 
should be abandoned. Rather, it is necessary to specify 
more precisely how to time the quantum particle in the 
tunneling regime. Different specifications lead to differ- 
ent relevant time scales; in other words, the importance of 
the differently defined times rests on the context where 
they become physically significant quantities. For ex- 
ample, the traversal time of Biittiker and Landauer (BL 
time) J!]], r = L/v sc , given by the barrier length L divided 
by the "semiclassical" velocity v sc — [2(V — i?)/m] 1 / 2 , 
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marks the transition from sudden to adiabatic regimes 
for an oscillating barrier jl] , and determines the rotation 
of the spin in a weak magnetic field in opaque conditions 
||; whereas the average over wavepacket components of 
the (monochromatic) "phase times" provides the mean 
arrival time of the transmitted wave packet [M ^) . These 
two time scales may be very different. The Biittiker- 
Landauer time increases with decreasing energies up to a 
finite value, whereas T Ph (0, L) (the so called extrapolated 
phase time) diverges as E — > 0, and tends for increasing 
L to a constant value, 2Ti/[v sc {2mE) 1 / 2 ]. This later prop- 
erty implies that the arrival of the transmitted wave be- 
comes independent of L (Hartman effect Jul), although 
the independence only holds until a certain critical length 
L c |7) where above-the-barrier components start to dom- 
inate. For L > L c the mean arrival time depends on L 
linearly. While r and T Ph are surely the most frequently 
invoked tunneling times, they do not exhaust all timing 
questions. 

In particular, in Ref. |pd| , the time t tr = t tr (x) that 
characterizes the transition from the transient to the sta- 
tionary regime for a wave formed by a point source with 
a sharp onset, was identified and the time of arrival, 
tp = tp{x), of the peak of the forerunner at a point x 
in opaque conditions was also identified. Here "opaque" 
means xkq ^> I, where kq = [2m(V — E^)] 1 / 2 /%, and 
Eq = fwjQ. This particular time scale turned out to be 
(surprisingly) proportional to the BL time r i.e., 
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even though the time-frequency analysis of the forerun- 
ner showed that it was composed by frequencies above 
threshold. In other words, it corresponded to non- 
tunneling. In spite of its frequency content, the fore- 
runner's peak "travels" with a velocity proportional to 
v sc , v p = 3 1 / 2 w sc , which increases with decreasing ener- 
gies, and its intensity diminishes exponentially as it pro- 
gresses along the coordinate x. Other works had already 
pointed out the dominance of non-tunneling components 
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in the forerunner |T^, [TJ], [TJ], 15, [T(| but had not char- 
acterized its time dependence. Their main objective was 
to show that a previous prediction by Stevens |l7j] (who 
believed that a tunneling monochromatic o^o-front, asso- 
ciated with a pole at u>o in the complex frequency plane 
would arrive at r) did not hold because of the effect of the 
saddle point or other critical points, such as resonances 
in a square barrier model [pL6[| . (For the source- with-a- 
sharp-onset model the frequency of the saddle is given by 
uj s = E s /~h, where the saddle energy, E s = V + x 2 m/2t 2 , 
coincides with the energy of a classical particle travel- 
ing from the source to i in a time t.) Biittiker and 
Thomas proposed to enhance the importance of the 
monochromatic front associated with loq compared to the 
forerunners by limiting the frequency band of the source 
or of the detector; it was shown later in |nj that the 
monochromatic front could not be seen in opaque condi- 
tions even with the frequency band limitation. However, 
a clear separation of the amplitude into two terms, one 
associated with saddle and forerunner and the other with 
the pole and the "monochromatic front" , is only possible 
for opaque conditions. Non semiclassical conditions have 
been much less investigated []l9|, even though these are 
actually easier to observe because of the stronger signal. 
In this paper we shall be mainly concerned with them. 

The motivation of the present paper is a recent publi- 
cation of two of us where the time evolution of an initial 
cutoff wave truncated at the left edge of a square barrier 
(shutter problem) f2(i[| was examined. There it was found 
that the probability density at the barrier edge x = L, 
exhibits at short times a transient structure named time 
domain resonance. The maximum, t p , of the time do- 
main resonance showed a plateau region for small L, or 
more accurately a shallow basin, followed by a linear de- 
pendence for larger L; this behavior is reminiscent of the 
Hartman effect, but the time of the plateau did not co- 
incide with the phase-time estimate. It was found that 
for a broad range of parameters, t p in the basin may be 
written approximately as, 
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where e\ and Eq, correspond to the energy of the first top- 
barrier resonance and the incidence energy, respectively. 
Such a dependence had not been described before. On 
the other hand, along the linear regime, at larger values 
of L, t p is described by, 
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(3) 



where V\ = hai/m, with a\ the real part of the first top- 
barrier pole k\ = a\ — ib\. The above time scales are 
also different from the BL-time, even though they may 
coincide with t^ for a particular value of L. 

We shall show in this paper that this basin dependence, 
and the corresponding time scale, is not only present at 
the transmission edge of the square barrier studied in p0[ . 



It may also be found, mutatis mutandis, in the sharp on- 
set source model examined in [[ll) for small xkq . We also 
find the basin for the time of the forerunner versus po- 
sition in the internal region of the square barrier, and 
for a step potential barrier. There are not resonances in 
the sharp onset source model, or for the step potential, 
so the time scale of the basin minimum is in these cases 
inversely proportional to Kq, namely to the difference be- 
tween the "potential level" V and the source main energy 
E . This is to be contrasted with the dependence on k^ 1 
of the traversal time r. 

An important open question was to determine if the 
forerunner at small lengths in the evanescent region cor- 
responds or not to tunneling frequencies. We shall show 
by a simple time-frequency analysis that the peak of the 
forerunner in this regime is composed predominantly by 
under-the-barrier components, so that indeed a genuine 
tunneling time scale different from phase or BL times has 
been found. 

In section II the source with a sharp onset is discussed, 
and in section III we shall examine the internal region of 
the square barrier, as well as the step barrier case. As 
the reader will soon discover, apart from certain peculiar- 
ities, all models show a small length region of the order of 
the penetration length k^" 1 where the forerunner is dom- 
inated by tunneling components and arrives at a time 
proportional to Kq 2 . 



II. SOURCE WITH A SHARP ONSET 

A. Formalism 

The time dependent solution of the Schrodinger equa- 
tion for a potential V(x) — V that occupies the entire 
space, subject to a source boundary condition, 

V>(0,t) = e- tuJot/h O(t), (4) 

has the form, |ll| 

t) = I e -W+« 2 /(4C* t ) [w{ _ u , o ) + w{ _ u n )} (5) 



where w(z) — e z erfc(— iz) and 
1 + i 
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t = x/v sc = xmj Koh, 
C = (H/2m) 1/2 . 



(6) 



(7) 



In the opaque limit (xkq 3> 1), an excellent approxima- 
tion is given by adding up the pole and saddle terms, 



$(x,t) pa *o+* s , 

* = e- luJat e- KaX e{t -t), 
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FIG. 1: |*| 2 (solid line), |* s | 2 (circles), |* | 2 (dashed line), 
and |*o + * s | 2 (dotted-dashed line) versus t for E = 0.9071/, 
and x = 2.75 nm. Here and in Figs. §|J V = 0.3 eV. In all 
figures m = 0.067m e , where m e is the electron's mass. 



For opaque conditions at fixed x the saddle term is totally 
dominant in the forerunner peak; at later times it fades 
away and the pole (tunneling) term becomes dominant. 
This allowed to identify t p from the peak of | 2 , and t tr 
as the time where the two contributions were of equal im- 
portance. The small xkq regime did not permit the same 
approximation treatment, so that t p could not be char- 
acterized by the formula found for the opaque case, Eq. 
(|l|). It may be argued however, that the small xkq case 
is in fact more interesting since it allows non-negligible 
signals (densities), whereas for large opacity the signal 
becomes exponentially small and very hard to detect. 



B. Examples 



FIG. 2: ujav/^v (solid line), luo/luv (dashed line), and uj 3 /ujv 
(circles) versus t for the same parameters of Fig. |l[ The 
square marks the value of t p . 
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In Fig. |l| we have depicted the density versus t for a 
fixed x smaller than the penetration length l//vo for the 
source-with-a-sharp-onset model. In all calculations the 
effective mass of the electron is taken as m = 0.067m e . 
The forerunner is identified as a smooth broad bump with 
its peak at t p . The contributions of the saddle and pole 
terms are also drawn, as well as the contribution of their 
sum. The saddle term reproduces the density only for 
very short times while the pole gives the large time be- 
havior. 

In Fig. |^ the average local frequency |nj ^l), ui av = 
— Im [(d^f /dt)/^f], the frequency of the saddle u> s , and ojq 
are shown, relative to wy, versus t for the same value of 
x. Note that uj av tends to u>o for large t, and to u) s for 
very short times, going from oj s to luq during the time 
span of the forerunner bump (compare with the den- 
sity in Fig. [I]); in particular uj av /ujy < 1 at the peak 
t p . Fig. @ shows the dependence with x of the time of 
arrival of the peak of the forerunner, t p , for two differ- 



FIG. 3: time of arrival of the forerunner's peak versus x for 
cjo/ = 0.907<x>v (upper solid line), and ojq = 0.544tjy (lower 
solid line). The upper and lower dashed lines are the corre- 
sponding BL-times. The values of 1/fto are also shown with 
a square {ujq — 0.544ojv) and a circle (luq — 0.907luv) 



ent values of ujq. The (unsharp) transition between the 
quasi-plateau region at lower x and the linear regime for 
larger x occurs around a few penetration lengths 1/ko- 
Thus, a smaller ujq implies a more reduced plateau. Also 
shown are the corresponding Biittiker-Landauer times for 
comparison. In the plateau region the forerunner's peak 
does not travel but arrives roughly simultaneously at all 
x (in fact up to the basin minimum the forerunner ar- 
rives earlier at larger a;!); beyond the plateau region the 
forerunner, whose peak is characterized in that case by 
Eq. ([j]), arrives later at larger x, and becomes dominated 
by frequencies "above the barrier" . 
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FIG. 4: uiav/ujv (solid line), ujo/lov (short dashed line), and 
uj s /ijjv (circles) for the same cases of Fig. ||. Upper curves: 
ojq — 0.907uiv', lower curves: ujq = 0.544o;y. The values of 
are indicated as in Fig. H 




FIG. 5: Density |*| 2 versus x for lo /lo v = 0.907 and t = 
0.05, 24, 48, 72 and 97 ns (solid, short-dashed, long-dashed, 
dot-dashed, and circles respectively). 



Fig. H shows the average frequency uj av calculated right 
at the forerunner's peak at the time t p (x) for each value 
of £, as well as u s , and ujq, all relative to cuy, for the 
same two values of uj of the previous figure. The fore- 
runner's peak frequency begins at ujq at x — > 0, then 
it grows smoothly and becomes non-tunneling (u> av /u>v 
crosses 1) around 1/kq. For larger x, u> av tends to the 
saddle frequency u s . Further insight into the nature of 
the tunneling forerunner is obtained by looking at a sc- 
ries of snapshots of the wave density versus x (instead of 
looking at density versus time for a fixed location). The 
plateau, or rather shallow basin, must be, according to 
Fig. (^) , the result of a "breathing" transient mode of the 
evanescent wave, which implies an essentially simultane- 
ous growth of the wave for all i up to a few penetration 
lengths, and until the maximum is achieved at t p . This is 
indeed confirmed in Fig ||, where several snapshots of the 
density are taken at five different times. In the final one 
the wave is already close to its asymptotic (large time) 
form. The linear dependence of the time of arrival of the 
peak with Kq 2 is shown in Fig. ^, which is drawn by 
choosing the values of V — Eq and t p corresponding to 
the minima of the shallow basins obtained for different 
values of E (see Fig. ||). 



III. QUANTUM SHUTTER 

A. Formalism 

Assume now an arbitrary potential V(x) (0 < x < L) 
that vanishes outside the internal region, and the initial 
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FIG. 6: Time of arrival of the forerunner's peak t p versus 
1/(V — -Bo) at the position that minimizes t p for a given Eo, 
see Fig. H. As in all the figures of this section, V = 0.3 eV. 



condition, 

*<,M = 0> = {^- : -'*-- (12) 

which corresponds to a plane wave impinging on a per- 
fectly reflecting shutter placed at x = 0, just at the left 
edge of the structure. The time dependent process begins 
with the instantaneous opening of the shutter at t = 0, 
enabling the incoming wave to interact with the poten- 
tial at t > 0. The exact solution along the internal region 
(0 < x < L) is given by ||, 

OO 

* l = <t> k M{y k ) - <t>- k M{y- k ) - ^ p n M(y kn ). (13) 

Tl— — OO 
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FIG. 7: (a) Snapshots of |* l | 2 (solid line) taken at different 
values of time: t\ = 1.0 fs, ti = 2.0, and t% — 4.0 fs. The 
stationary solution \4>k\ 2 is included for comparison, and a full 
square indicates the value of Kq , which is the penetration 
length of the stationary solution, (b) Time evolution of |\l/ l | 2 
at different fixed values of the position x: x\ — 0.5 nm, X2 = 
0.7 nm, X3 — 1.0 nm, and x^ — 2.0 nm. The solution exhibits 
a behavior similar to diffraction in time pattern. 



In the above expressions, the quantities <f>±k refer 
to the stationary solution, and the factor p n = 
2iku n (0)u n (x) I (k 2 — fc 2 ) is given in terms of the reso- 
nant eigenfunctions, {u n (x)}, with complex eigenvalues 
p2[ k n = a n — ib n (a n ,b n > 0). The index n runs over 
the complex poles k n , distributed in the third and fourth 
quadrants in the complex fc-plane. Both the complex 
poles {kn} and the corresponding resonant eigenfunc- 
tions {u n (x)}, can be calculated using a well established 
method, as discussed elsewhere |^0|, p2| . In the above 
equation the M functions are defined as, 

M{y q ) = ^w(iy q ), (14) 

where y q = —e~ l7T / 4 [(h/2m)t] 1 / 2 q, q = ±k, and k± n . 



FIG. 8: Exact calculation of t p (solid dots) as a function of 
the position x. The parameters are the same as in Fig. [?]. A 
full square indicates the value of Kq 1 . 



B. Examples 

In what follows we shall explore the features of the 
probability density along the internal region of one- 
dimensional rectangular potential barriers of height V 
and thickness L, defined along the interval < x < L. 
We choose the following parameters: V — 1.0 eV, L = 
40.0 nm, and incidence energy E$ = h 2 k 2 /2m = 0.1 
eV. For this particular case the full barrier is opaque, 
a = LftmV] 1 / 2 /h = 53.0, but we are interested in dis- 
tances of the order of, or smaller than, the penetration 
length Kg 1 << L. In Fig. (a) we use Eq. (|l||) to 
explore the birth of the probability density along the in- 
ternal region of the potential. Here we plot |\I> 1 | 2 (solid 
line) as a function of small values of the position x, for in- 
creasing values of time. Notice that as time increases, the 
probability density evanesces along the internal region 
swinging around the stationary solution (dashed line). 

In Fig. (b) we illustrate the time dependence of the 
transient "swinging mode" exhibited in Fig. [?] (a). We 
exhibit the behavior of l^l 2 (solid line) as a function of 
time for several values of the position xi (i — 1, 2, 3,4). 
The probability density grows monotonically until it 
reaches its first maximum value, and oscillates thereafter 
as it tends to the stationary situation given by |(/)fc| 2 [p2[. 
Here we observe again a peculiar behavior in the time 
evolution of l^*! 2 , presented in plot (b) of Fig. [?]: the 
first maximum of the probability density appears earlier 
at larger values of x. In Fig. |^ we plot the exact calcu- 
lation of t p (solids dots), which corresponds to the peak 
value of the first maximum of | \I/ * | 2 , measured at different 
values of the position, x. 

In Fig. |^ we can clearly identify two regimes, as in 
the case of Fig. ||. The first of them corresponds to a 
basin, where for small values of the position, t p decreases 
reaching a minimum value as x increases. This is in ac- 
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FIG. 10: Maximum of the time domain resonace, t„, as a 
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FIG. 9: Plot of the relative local frequency ujav/cov (solid 
line) for the positions (a) x\ — 0.5 nm, , and (b) X4 — 2.0 
nm, corresponding to the cases depicted in Fig. m (b). The 
cutoff-frequency luv /ofy = 1 (dotted line) is included for com- 
parison, and a full square indicates the position of t p . 



cordance with the observed behavior of the maximum in 
Fig. |?] (b). However, if x is further increased, t p begins to 
grow as a function of the position. Apparently this sec- 
ond regime corresponds to a situation where t p increases 
linearly with x. The result depicted in Fig. |^ appears 
to be a replica of a similar behavior of the probability 
density at the barrier edge x = L, recently reported in 
Rcf. 

We have also analyzed the frequency content of l^l 2 
along the internal region. In Fig. ^| (b) we plot u) av /u>v 
along the relevant time interval for some of the values of 
the position, depicted in Fig. ^ (b). We can appreciate 
that the probability density at this small values of the 
position x, is composed entirely by under-the-barrier fre- 
quency components i.e oj av /uJv < 1 in the vicinity of the 
maximum t p . Notice that this result is similar to the one 
obtained in Fig. p| 

Finally, Fig. 10 shows the linear dependence of the 



function of 1/(V - E ). 
and L = 40.0 nm. 



The parameters are the V = 1.0 eV, 




FIG. 11: LO av /iOv (solid line), uq/ljv (dotted line), and lu s /lov 
(dashed line) for V = 1 eV and Eg = 0.5 eV. The value of 
Kg 1 is indicated with a square. 



l/(V — E ) in analogy to the analysis of the point source 
problem. Note, however, that in our case we could also 
have considered a plot with respect to l/(e\— E ) with a 
similar result since it turns out that for opaque barriers 
£1«V. 



C. Potential step 

The potential step VQ(x) requires a different formal 
treatment |2l| but very similar results are found. We 
have evaluated (x\^>(t)) by means of an integral over the 
energy eigcnfunctions for the the initial state, 



time corresponding to the basin minimum with respect to 



^(x,k:t = 0) = 



,,tkx 



-OQ < X < 0, 

x > 0. 



(15) 
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FIG. 12: Time of arrival of the forerunner's peak t v versus 
1/(V — Eo) at the position that minimizes t p for a given Eo- 
V = 1 eV. 

In Fig. [ll] the average frequency peak of the forerun- 
ner is shown versus x for this system. As in the source 
with a sharp onset case, it goes from oj to oj s , cross- 

fthe threshold frequency at approximately 2/kq. Fig. 
shows also the linear dependence of the time of the 
forerunner (calculated at its minimum value for a given 
energy, i.e., at the basin minimum) versus 1/(V — Eq). 

IV. FINAL COMMENTS 

By using different time-dependent models of quantum 
waves impinging on a potential barrier, we have shown 



the existence of forerunners dominated by under-the- 
barrier components. This transient structures are indeed 
associated with a true tunneling process, occurring at 
distances smaller than the penetration length, Kq 1 . We 
find that the time of arrival is proportional to 1/(V — Eo), 
where Eq stands for the incidence energy, chosen below 
the barrier height, V. This is in contrast to the dom- 
inance of above the barrier frequencies for larger dis- 
tances, and to other well known tunneling time scales 
such as the Biittiker-Landauer or "phase" times. 

These results suggest a number of questions that need 
further investigation. For example, the proportionality 
constant observed in the arrival times of the forerunner's 
peak, depends on the peculiarities of each model, and a 
general theory for its specific value should be found. Also, 
the Larmor time for spin rotation defined by Biittiker is 
given for very thin barriers by h/V in the limit k — [^). 
In spite of the very different ways in which the Larmor 
time and the forerunner's peak time are obtained, these 
two quantities might be related. Finally, we hope that the 
results presented may motivate an experimental search of 
tunneling forerunners and their corresponding time scale. 
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